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> ■ ABSTRACT 

in 

CO . We report the discovery of water vapour toward the carbon star V Cygni. We have used Herschel's HIFI instrument, in dual beam 

04 ' switch mode, to observe the In - Ooo para-water transition at 1113.3430 GHz in the upper sideband of the Band 4b receiver. The 

observed spectral line profile is nearly parabolic, but with a slight asymmetry associated with blueshifted absorption, and the integrated 
antenna temperature is 1.69 + 0.17 K km s"'. This detection of thermal water vapour emission, carried out as part of a small survey 
of water in carbon-rich stars, is only the second such detection toward a carbon-rich AGB star, the first having been obtained by the 
5^ ' Submillimeter Wave Astronomy Satellite toward IRC+I0216. For an assumed ortho-to-para ratio of 3 for water, the observed line 

intensity implies a water outflow rate ~ 3 - 6 x 10"^ Earth masses per year and a water abundance relative to H2 of ~ 2 - 5 x 10"* . This 
value is a factor of at least lO'* larger than the expected photospheric abundance in a carbon-rich environment, and - as in IRC+ 10216 - 
raises the intriguing possibility that the observed water is produced by the vapourisation of orbiting comets or dwarf planets. However, 
. J ■ observations of the single line observed to date do not permit us to place strong constraints upon the spatial distribution or origin of 

Itn] ' the observed water, but future observations of additional transitions will allow us to determine the inner radius of the H20-emitting 

M zone, and the H2O ortho-to-para ratio, and thereby to place important constraints upon the origin of the observed water emission. 

C^ 

Key words. Stars: AGB and post- AGB - Stars: circumstellar matter - Submillimeter: stars 

1 . Introduction in photospheric composition, water vapour has previously been 

discovered, using the Submillimeter Wave Astronomy Satellite 

The carbon-to-oxygenratio is the critical determinant of the pho- (SWAS), in the circumstellar outflow of the extreme carbon star 

tospheric chemistry m evolved stars; the photospheres of oxygen irc+10216 (Melnick et al. 2001), for which the C/O ratio is 

rich-stars, with C/O ratios < 1, are dominated by CO and H2O, ^ 14 xhe observation of water vapour in fliat source, with 

while those of carbon-rich stars are dominated by CO, HCN, and ^^ inferred abundance relative to H2 of ~ 10"^ (Agundez & 

C2H2 and contain very litde H2O. Despite this sharp dichotomy Cernicharo 2006; Gonzalez- Alfonso, Neufeld & Melnick 2007), 

has led to the suggestion of several possible origins for the water 

* Herschel is an ESA space observatory with science instruments pro- vapour, including (1) the vapourisation of icy objects (comets 

vided by European-led Principal Investigator consortia and with impor- or dwarf planets) in orbit around the star (Ford & Neufeld 
tant participation from NASA. 
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2001); (2) Fischer-Tropsch catalysis (Willacy 2004); (3) pho- 
tochemistry within an outer, photodissociated shell (Agundez 
& Cernicharo 2006); (4) photochemistry within a clumpy out- 
flow (Decin et al. 2010). In both carbon- and oxygen-rich stars, 
Cherchneff (2006) and Decin et al. (2008) have argued for the 
importance of shocks in setting non-equilibrium abundances of 
the molecules CO, HCN, CS, and SiO within the inner envelope, 
although existing models do not appear to predict a persistent 
enhancement of H2O in stars with C/O > 1 . 

Using the HIFl instrument (de Graauw et al. 2010) on 
Herschel (Pilbratt et al. 2010), it is now possible to search with 
unprecedented sensitivity for water vapour in carbon-rich stars. 
As part of the HIFISTARS Guaranteed Time Key program, we 
will carry out a survey for water vapour in eight carbon-rich stars 
at distances greater than that (~ 170 pc) of IRC -(-10216. Here, 
we report the detection of water vapour in the very first carbon 
star targeted in this survey: V Cygni, a Mira variable of spectral 
type C6 (Wallerstein & Knapp 1998), with an apparent V mag- 
nitude varying from 13.9 to 7.7 mag, a period of 421 days, and a 
distance of ~ 400 pc (Bieging & Wilson 2001). 

2. Observations and results 

We observed the lowest transition of para-water, the In - Ooo 
transition at 1113.3430 GHz line, in the upper sideband of the 
Band 4b HlFl receiver. This observation, of duration 1 167 s in- 
cluding overheads, was carried out on 2010 March 5, using the 
dual beam switch (DBS) mode and the Wide Band Spectrometer 
(WBS). The WBS has a spectral resolution of 1.1 MHz, corre- 
sponding to a velocity resolution of 0.30 km s"^ at the frequency 
of the 1 1 1 - Ooo transition. The telescope beam, of half-power- 
beam-width (HPBW) ~ 20", was centered on V Cygni at coordi- 
nates a = 20h41m 18.3s, 5 = 48° 08' 29"(J2000). The reference 
positions for this observation were located at offsets of 3' on ei- 
ther side of the source. 

The data were processed using the standard HlFl pipeline 
to Level 2, providing fully calibrated spectra of the source. The 
Level 2 data were analysed further using the Herschel Interactive 
Processing Environment (HIPE; Ott 2010), version 2.4, along 
with ancillary IDL routines that we have developed. Having 
found the signals measured in the two orthogonal polarizations 
to be in excellent agreement, we combined them to obtain an av- 
erage spectrum. Figure 1 shows the WBS spectrum of para-H20 
1 11 - Ooo obtained toward V Cygni, with the frequency scale ex- 
pressed as Doppler velocities relative to the Local Standard of 
Rest (LSR) and the intensity scale expressed as antenna temper- 
ature. A zeroth-order baseline has been subtracted. The vertical 
dashed line indicates the LSR velocity of the source, as deter- 
mined by Bieging & Wilson (2001, hereafter BWOl) from ob- 
servations of the CO 7 = 2-1 line. The integrated antenna tem- 
perature is found to be 1 .69 ± 0. 17 K km s"' . Our identification 
of the observed feature with the para-H20 1 1 1 - Ooo transition is 
supported by the absence of any plausible alternative candidate 
in the JPL (Pickett et al. 1998) or CDMS (Muller et al. 2001) 
spectral line catalogues. 

3. Derivation of the water outflow rate and 
abundance 

In modeling the water emission from V Cygni, we have used the 
methods described by Gonzalez-Alfonso, Neufeld & Melnick 
(2007, hereafter GNM). To determine the H2O rotational pop- 
ulations and the resultant H2O emission spectrum, we included 
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Fig. 1. Continuum-subtracted spectrum of para-H20 In - Ooo 
obtained toward V Cygni. The vertical dashed line indicates the 
LSR velocity of the source, as determined by BWOl from obser- 
vations of the CO 7 = 2-1 line. 

Table 1. Assumed parameters for V Cygni 



Distance: 400 pc 
Stellar luminosity: 8500 Lq 
Effective temperature: 2500 K 
Stellar radius: 3.42 x 10" cm 
Systemic LSR velocity: -1-13.8 km s"' 
Terminal outflow velocity: 1 1.8 km s"' 



the effects of radiative pumping - through the V2 = 1 and V3 = 1 
vibrational states - and of collisional excitation by H2, together 
with a treatment of radiative transfer based upon that discussed 
by Gonzalez-Alfonso & Cernicharo (1997). 

The assumed parameters for the source are summarized in 
Table 1 and discussed below. A variety of distance estimates for 
V Cygni have appeared in the literature. BWOl obtained values 
of 406 pc and 456 pc respectively by using the P - K (period 
- absolute K magnitude) and P - L (Period - bolometric lumi- 
nosity) relationships presented by Groenewegen & Whitelock 
(1996; hereafter GW96). However, our own fit to the spectral 
energy distribution (SED), shown in Figure 2, implies that a dis- 
tance of 342 pc is required to fit the GW96 P - L relationship. 
Here, we used a dust radiative transfer model to fit a combination 
of flux measurements from Tycho, Hipparcos, MSX, 2MASS, 
ISO, IRAS, and JCMT. This model, full details of which will be 
presented in a future publication (Schmidt et al. 2010), assumes 
a power-law distribution of grain radii with index -3.5 between 
0.005 and 0.25 /im, a maximum dust temperature of 1050 K, an 
inner radius for the dust shell of 2 x lO'^ cm, and a total opti- 
cal depth Ay = 9.45 mag. A smaller distance estimate (370 pc) 
than those of BWOl was also obtained by Schoier & Olofsson 
(2000). All of the values mentioned above are consistent with 
the Hipparcos-measuied parallax of 3.69 + 1 .77 mas, which im- 
plies a 1 cr distance range of 180 - 520 pc. Accordingly, we adopt 
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Fig. 2. Spectral energy distribution of V Cygni. A dust radiative 
transfer model (see text) has been used to fit the fluxes measured 
by Tycho, MSX, 2MASS, ISO, IRAS, and JCMT. The dashed 
curve shows the unattenuated photospheric emission, while the 
solid curve shows the total emission from the attenuated photo- 
sphere and surrounding dust. 



a value of 400 pc for the estimated distance of V Cygni, with a 
likely uncertainty ~ 20%. 

This estimate of the distance then requires a bolometric lu- 
minosity of 8500 Lq to match the observed SED shown in Figure 
2, and a stellar radius of 3.42 x 10'^ cm, given an assumed ef- 
fective temperature of 2500 K. The values given in Table 1 for 
the systemic velocity of the source and the terminal velocity of 
the outflow were based upon Herschel/HTPl observations of the 
CO y = 6 - 5 ti-ansition at 691.473 GHz (Schmidt et al. 2010) 
and are in excellent agreement with previous determinations by 
BWOl. To obtain an estimate of the total mass-loss rate in the 
outflow, we have modeled the fluxes observed by Herschel/HIFI 
fortheCOy = 6-5, 7= 10 -9 and 7= 16- 15 transitions. Full 
details of the observations and modeling of CO will be given by 
Schmidt et al. (2010). Our best estimate of the gas mass-loss rate 
in the inner envelope is 4.6 x lO^^Moyr"', a factor ~ 3 larger 
than that obtained by Schoier & Olofsson (2000) from a fit to 
the CO 7 = 2-1 transition, and the derived CO/H2 ratio is 10"^. 
Indeed, we find that a constant mass-loss rate model that fits the 
CO 7 = 6-5, 7= 10-9 and 7=16-15 transitions substan- 
tially overpredicts the flux in the CO 7 = 2-1 transition. This 
discrepancy suggests some variability in the mass-loss rate, with 
a larger value applying to the inner envelope where the higher- 
lying transitions of CO originate. We obtained a satisfactory fit 
to both the CO rotational line fluxes and the continuum spectrum 
by assuming a gas and dust density that decreases as radius"^'^, 
instead of the radius"^ density profile expected for an envelope 
with a constant mass loss rate and outflow velocity. The gas-to- 
dust mass ratio in this model is 510. 

Fortunately, given the significant uncertainties in many of the 
assumed parameters listed in Table 1, the derived water outflow 
rate is not strongly dependent upon any of them. As discussed 
in GNM, unless the mass-loss rate is extremely large, the exci- 
tation of water is dominated by radiative pumping via the 6 fim 
V2 band. Thus, for a given water outflow rate, the observed water 
line fluxes scale linearly with the observed 6 //m continuum flux. 
Since the latter is an observed (rather than a derived) quantity. 
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Fig. 3. Continuum-subtracted and smoothed spectrum of the 
1 1 1 -Ooo transition of para-water at 1 1 13.3430 GHz (black), com- 
pared with the the best fit from Model B (blue). The residuals are 
shown in red. Here, the vertical axis is the main beam brigthness 
temperature, obtained by dividing the antenna temperature by an 
assumed main beam efficiency of 0.67 



our estimate of the water outflow rate is largely independent of 
the distance or total outflow rate assumed for the source. 

In modeling the water line strength and profile observed to- 
ward V Cygni, we have investigated two models for the spatial 
distribution of the observed water In Model A, we assume that 
water is present at radii as small as 4.5 x lO''* cm, while in Model 
B, we adopt an inner radius Ri„ =2x10'^ cm at which the water 
is injected into the outflow. In both cases, we assume an outer 
radius of /?out = 1 x 10'* cm, the estimated photodissociation 
radius for waterQ 

The distribution assumed in Model B is expected if the 
vapourisation of icy objects is the origin of the observed wa- 
ter vapour, since all such objects at smaller distances from 
the star will have been vapourised already (Ford & Neufeld 
2001, GNM). In both models, we adopted the profile given by 
Goldreich & Scoville (1976) for the outflow velocity as a func- 
tion of radius. For Model B, this velocity profile was found to 
give a better fit to the observed line shape than an alternative 
velocity profile, based upon the work of Skinner et al. (1999), 
which entails a more rapid acceleration of the outflowing mate- 
rial. A clear asymmetry in the line profile, apparently associated 



' The mean lifetime against photodissociation for a water molecule 
in the unshielded mean Galactic ultraviolet radiation field is 1.2 x 10' s 
(Roberge et al. 1991), corresponding to a distance of 1.4 x 10'^ cm at 
the terminal outflow velocity of the V Cygni outflow. Shielding by dust 
in the outflow increases this distance by a factor of several, leading to 
photodissociation radii ranging from ~ 5 - 10 x 10'' cm"^ for radiation 
fields within a factor of two of the mean Galactic value. We adopt the 
upper end of this range, in order to obtain a conservative estimate of 
the water outflow rate. To account for the observed water line strength, 
smaller values of R„„t would require larger values of the water abun- 
dance. 
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Table 2. Derived parameters for V Cygni 



Table 3. Predicted line strengths for V Cygni 



Model A: Water inner radius = 4.5 x 10" cm 

Para-water abundance relative to H2: 5.3 x 10"^ 
Total water abundance" relative to H2: 2.1 x 10"'' 
Total water outflow rate: 3 x 10"' Earth mass / yr 

Model B: Water inner radius = 2 x 10" cm 

Para- water abundance relative to H2: 1.0 x 10"* 
Total water abundance" relative to H2: 4.2 x 10"* 
Total water outflow rate: 6 x lO""" Earth mass / yr 
"for an assumed ortho-to-para ratio of 3 



with blueshifted absorption, suggests a significant turbulent ve- 
locity width, of order 1.3kms "'. Equally good fits can be ob- 
tained for models A and B, implying that our current observa- 
tions (of a single line) cannot be used to provide useful con- 
straints on the spatial distribution or origin of the water vapour 
in V Cygni. The best-fit parameters are given in Table 2, for each 
of the two models, and the fit for Model B is shown in Figure 3. 
For an assumed water ortho-to-para ratio of 3, the required water 
outflow rates are 3 x 10"^ and 6 x 10"^ Earth masses per year for 
models A and B respectively. 



4. Discussion 

Our detection of thermal water vapour from V Cygni is only 
the second such detection from the circumstellar envelope of 
a carbon-rich AGB star, the first having been obtained toward 
IRC+10216. The water outflow rate in V Cygni, 3 - 6 x 10"^ 
Earth masses per year, is 4 - 8 times that in IRC+10216, even 
though the total mass outflow rate of 4. 6x10"'' Mq yr"' is smaller 
by a factor ~ 8. Thus the water line luminosity in V Cygni im- 
plies a water abundance ~ 30-60 times as large as that inferred 
for IRC-H10216. The large H2O abundance -3-6x10"'' that 
we infer is at least a factor of 10"* larger than that predicted in 
models for the photospheres of C-rich stars in thermochemi- 
cal equilibrium (e.g. Cherchneff 2006). As in IRC+10216, the 
inferred abundance raises the intriguing possibility that the wa- 
ter is produced by the vapourisation of orbiting comets or dwarf 
planets. In this picture (Ford & Neufeld 2001), a collection of 
orbiting icy objects in a Kuiper-belt analogue are vapourised as 
the star ascends the asymptotic giant branch and its luminosity 
gradually increases; to match the observed water outflow rate in 
V Cygni, this scenario would require a reservoir of ~ 10 Earth 
masses of water ice available for release during the AGB phase. 
Observations of the single water line observed to date, how- 
ever, do not permit us to place strong constraints upon the spatial 
distribution or origin of the observed water, but future observa- 
tions of additional transitions will allow us to determine the inner 
radius, R^a, of the HaO-emitting zone, as well as the H2O ortho- 
to-para ratio. The radius, /?„, is a critical diagnostic, because 
the comet vapourisation model predicts a value ~ 2 x lO'^cm 
(Model B), all orbiting icy objects within that radius having al- 
ready been vapourised. In Table 3, we present predictions for 
the strengths of various emission lines lying within the spectral 
range accessible to the HIFI, SPIRE and PACS instruments on 
Herschel. In the HIFI (and SPIRE) spectral range (upper section 
of Table 3), results are given as integrated antenna temperatures; 
in the PACS spectral range (lower section), the values given are 
fluxes, in units 10"'*" W m"^. The predictions presented here are 



Transition 


Frequency 


Wavelength 


Line strength 




(GHz) 


Oum) 


/ Tmb dv i 
Model A 


nKkms"' " 
Model B 


1 1 1 - Ooo 


1113.343 


269.272 


2.53" 


2.53* 


2o2 - In 


987.927 


303.456 


2.46 


2.79 


2ii - 2o2 


752.033 


398.643 


1.18 


1.17 


220 - 2ii 


1228.789 


243.974 


0.56 


0.56 


422 - 331 


916.171 


327.223 


0.47 


0.02 


lio - loi 


556.936 


538.289 


1.16 


1.06 


2l2 - loi 


1669.905 


179.527 


6.40 


5.90 


221 - 2l2 


1661.008 


180.488 


1.92 


1.79 


3o3 - 2l2 


1716.770 


174.626 


5.43 


4.93 


3l2-221 


1153.126 


259.982 


1.37 


0.85 


3 12 - 3o3 


1097.365 


273.193 


1.67 


1.49 


321-3,2 


1162.911 


257.795 


0.65 


0.59 


432 - 5o5 


1713.882 


174.920 


0.64 


0.03 


Transition 


Frequency 


Wavelength 


Line strength 




(GHz) 


(pm) 


(10"'*Wm"2) 








Model A 


Model B 


2ii-2o2 


2968.748 


100.983 


1.27 


0.82 


331 - 220 


4468.572 


67.089 


1.16 


0.42 


221 - lio 


2773.976 


108.073 


2.20 


1.67 


321-212 


3977.045 


75.381 


2.61 


1.11 


4l4 - 3o3 


2640.473 


113.537 


1.70 


0.70 


330 - 221 


4512.385 


66.438 


2.51 


1.12 


432 - 3l2 


3807.256 


78.742 


1.72 


0.23 


5o5 - 4l4 


3013.199 


99.493 


1.23 


0.22 


432 - 321 


5107.284 


58.699 


1.73 


0.35 


616 - 5o5 


3654.602 


82.031 


1.20 


0.02 



"integrated main beam temperature 
''observed line strength 



for an assumed ortho-to-para ratio of 3. Several lines, such as 
the 916.171 GHz 422 - 33i transition observable by HIFI and 
the 616 - 5o5 82.031 /im transition observable by PACS, show 
a strong dependence upon Ria and should provide an excellent 
discriminant between difi'erent models. 
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